Probing dark matter and CMSSM with same-sign dilepton searches at the LHC 
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We introduce new observables for the study of the inclusive same sign dileptons production at 
LHC which are built out of ratios of the observed number of same-sign dileptons, both with same 
N(£,£) and different flavor N (£,£'). As a case study we apply them to the stau coannihilation 
region of the constrained minimal supersymmetric standard model. We show that the new variables 
depend rather mildly on the center of mass energy and how these can be used to constraint the 
parameter space in the (mj/j, tan /3) plane. 
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Introduction. The starting of the LHC era will al- 
low us to finally shed light into the last missing piece of 
the standard model (SM), the Higgs boson, and hopefully 
to probe the supersymmetric (SUSY) extension of the SM 
that represents the most popular solution to the hierar- 
chy problem, gauge coupling unification and the nature 
of dark matter. In the minimal supersymmetric standard 
model (MSSM) with R-parity conservation and in the 
mSUGRA inspired constrained version (CMSSM), the 
lightest neutralino is the lightest SUSY particle (LSP), 
neutral and stable and one of the favored dark matter 
(DM) candidates 0. 

In the framework of the SM, events in proton-proton 
(pp) collisions with two isolated same-sign leptons in the 
final state, or same-sign dileptons (SSD), are very rare. 
They may come from double gauge boson production 
WZ, WW and decays, double parton scattering or ttW, 
the last processes yet to be observed in proton-proton 
collisions. This makes this signature very natural to look 
for new physics. Isolated SSD represent also a standard 
search channel for SUSY models |J . The main processes 
that lead to inclusive final states SSD in proton-proton 
collisions are gluino pair production pp — > gg, gluino- 
squark associate production pp — > gq and squark pair 
production pp — > qq' . Cascade decays from these pairs 
produce easily XiXii X^X^ X2X2 that decay and even- 
tually lead to SSD that can be of the same flavor (SF- 



SSD) e ± e ± , 



e ± fJ , ± 



fX ± fX ± , 

e ± r ± 



r ± r ± or of different flavor (DFSSD), 
. The signal was searched for by CDF 
at Tevatron [| and both the CMS and ATLAS 0, i] 
collaborations have already performed searches for these 
particular class of events in the data sample of the final 
2011 run with the LHC working at the center of mass 
energy yfs = 7 TeV and a total integrated luminosity of 
about 5 fb _1 . No evidence for new physics was found 
and upper bounds on the number of SSD were used to 
set constraints on the parameter space of the CMSSM. 

In this brief report we show that by simply counting 
the number of SSD pairs N(ee), N(mj), N(tt), N(e/J,), 
N(e,r), N([it) it is possible to obtain direct information 
on fundamental CMSSM parameters like, for example, 
and tan/3. The proposed variables, that are build 



from the number of SSD, may give preliminary informa- 
tions on the fundamental parameters of the underlying 
theory. Once these are known, one has a guide to which 
decay chain are likely to show up in the data yields and 
thus measure the masses of the SUSY particles. As a 
practical example we study the SSD signature in connec- 
tion with the so-called stau co-annihilation region (tcr) 
of the CMSSM parameter space that is of interest for 
dark matter searches. 

Observables for the SSD signal. Let us consider 
gluino pair production. The gluino is a Majorana par- 
ticle and decays with equal branching fractions (B) to 
particles and antiparticles; this property allows to have 
two same sign charginos from the decay chains of the 
pair of sparticlcs produced in the pp collisions. The two 
same-sign charginos lead to final states with SSD, neu- 
trinos and LSP. Under this assumption, we introduce the 
new observables as follows. The cross section can be ap- 
proximated by cr(pp — > 2\l + X) ~ <r(pp —> gg) x B(g — > 
qq) 2 x B(q — > q'xt) 2 - We synthetically call this cross sec- 
tion a^ x and the branching ratios corresponding to the 
various chargino's decay chains leading to a lepton plus 
undetected particles (neutrinos and LSP) plus hadronic 
jets, xt ^ + X with Bi t t, I = e, pi, r. For a given inte- 
grated luminosity, the number of SFSSD is estimated as 
N(M) oc cr|| x {J2i B he) 2 , while the number of DFSSD 
is instead N{W) cx ^{ [Y.B,, ) x (Ei^O- Th c fac- 
tor 2 takes into account the fact that the leptons come 
from two identical charginos. In reason of the expected 
similar behavior of the first two lepton generations and 
the peculiar role held by the leptons and sleptons of the 
third family, we consider the ratios: 



N(It) 
N(tt) 

with I = e, n and 



2R, 



R = 



N(ll) 
N(tt) 



R 2 
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We remark that when all contributing mechanisms are 
considered the total cross section (<r xx 



a" 

xx 



a qg 

xx 



2 



. . . ) drops out exactly. Considering the relations ([!]) , 



2R, 
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(4) 
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we finally define the variables: 

1 JV(er) + iVQir) 
1 2 JV(tt) 

N2 = l MW) R ■ 

We choose the normalization factors 1/2 and 1/4 in such 
a way that JVi and N2 are average quantities that should 
take values similar to the basic ratios in Eq. (JXJ) . By 
definition, these variables should depend very mildly on 
the production cross section and on the center of mass 
energy of the machine. If we relax the assumption that 
the SSD are produced through a pair of charginos xtxt 
the above Eq. (j3]) and Eq. ((4]) are still valid to a very 
good degree of approximation although the cancellation 
of the cross section does not apply any longer. 

Dark matter and the CMSSM. The CMSSM is 
specified by assigning the value of the common gauginos 
mass mi/2, the common scalars mass mo, the common 
trilinear scalar coupling Aq at a certain unification energy 
scale. The other fundamental parameter is tan /3, the ra- 
tio of the two vacuum expectation values of the Higgs 
doublet. Further we assume the Higgs mixing parame- 
ter to be positive. The tqr is one of the regions Q 
of the parameter space in which the relic density of the 
neutralino as main dark matter component is compati- 
ble with the WMAP Q measurement of the cosmic mi- 
crowave background anisotropics, 0.096 < ilh 2 < 0.128 
at 3cr. In the tcr [9|, the neutralino relic density is con- 
trolled by co-annihilation before the freeze-out between 
the LSP and the lightest stau, f\, which is the next LSP 
with a mass splitting Am = — m Xo of few GeVs. 

In literature there exist theoretical studies discussing 
signatures of the tcr at LHC. In Refs. [l(| the authors 
show that with the reconstruction of the decays X2 ~> 
tit — > ttxi at the end of the gluinos and squarks cas- 
cade, it is in principle possible to measure the masses of 
the relevant particles and infer the values of the CMSSM 
parameters. In Ref. [ll[ it is argued that the measure- 
ment of the tau polarization in the fi — > txI can be use- 
ful to measure the mass difference Am. The smallness 
of Am is an important quantity for the phenomenology 
at LHC [l(|, as it sets the scale of the transverse mo- 
mentum of the leptons at the end of the cascade decays, 
but also for indirect detection signals in astrophysical 
searches with gamma ray, as it is responsible for spectral 
features at the endpoint of the gamma spectrum [12| and 
for the large cross section for the annihilation into a f\ 
pair at the galactic center (l3| . 

We consider the parameter space with fixed trilin- 
ear scalar coupling Aq = 0, which is used as bench- 
mark for supcrsymmetric studies at the LHC. The stau- 
coannihilation strips shown in Figure [T] are obtained 
with microOMEGAs [l4| using SoftSusy [r| as super- 
symmetric mass spectrum calculator, imposing WMAP 
constraints on the relic density accelerator constraints 
on the lightest Higgs, m^ > 114.4 GeV, chargino mass 
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Figure 1. Regions of the CMSSM parameter space in the 
(mojfB!^) plane which satisfy the WMAP constraint on the 
relic density 0.096 < Qh 2 < 0.128 for different values of tan/3 
that is given by the numbers above the strips. The dashed, 
dot-double-dashed and dot-dashed lines are respectively the 
the 95% CL exclusion curves obtained by CMS 0] and AT- 
LAS [1 H with the SSD search and the CMS SUSY search 
with hadronic final states 



m + > 103.5 GeV and the flavor physics constraint from 
bottom mesons decays B s —> A 7 and B s — > In 
Figure Q] we also report the 95% CL exclusion curves ob- 
tained by CMS and ATLAS HQ with the SSD search 
and the CMS SUSY search with hadronic final states [r| . 
As can be seen the first two years of operation of the LHC 
could only marginally exclude the parameter space of in- 
terest in this work. 

Analysis. To check the validity of the new observ- 
ables we make a preliminary simulation selecting points 
along the strips with tan j3 = 10 and tan (3 = 40. The the- 
oretical ratio R defined in Eq. @ is calculated in each 
point using SusyHIT [llj to compute the branching frac- 
tions. All informations of the selected models are then 
passed to Pythia 8.1 [l8[ to generate events in p-p colli- 
sions at s/s = 14 TeV. We have generated 4 x 10 6 events 
for each CMSSM point requiring that in the final state 
there are same-sign leptons, jets and missing energy. We 
find that Eqs. and (01 are well satisfied, the number 
of SSD coming from two same-sign xt ^ s correctly pre- 
dicted in terms of the theoretical ratio of the chargino's 
branching ratios R. When all production mechanisms 
are allowed, and especially at large tan/3, the number of 
taus SSD is contaminated by the decay chain involving 
the second neutralino, X1X2 an d X2X2 thus, depending 
on the point of the parameter space, deviations up to 
50% are observed. 

We now want to relate to TO1/2 and tan/3, the two 
CMSSM parameters that are the more interesting from 
the dark matter phenomenology point of view. In fact in 
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Figure 2. Contour maps in the plane (mi i%, tan 0) of the observables N\ (left column) and N2 (right column) as defined in 
respectively in Eq. Q and Eq. Q . Top plots are obtained with a simulation with yfs = 8 TeV, the bottom plots with y/s = 14 
TeV. The cuts and efficiencies are discussed in the text. 



Ref. |20J it was found that the neutralino mass along the 
strips of Fig. [I] is roughly given by m x = 0.44 x mi/2 — 
16 GeV for all the values of tan/3. Furthermore the spin- 
independent ncutralino-nucleon cross section, and hence 
direct detection rates, strongly depend on tan 0. 

In order to cover the full parameter space we took 20 
equally spaced points (in m 1 / 2 ) along each strip of Fig. Q] 
and the corresponding value of too is chosen to be the 
center value of the strip width along the mo direction. 
We carry out the simulation with 4xl0 6 events for each 
model point. For realistic results we employ the efficiency 
model for the SSD signal developed by the CMS collab- 
oration that allows to obtain realistic results at the gen- 
erator level bypassing the full detector simulation 0, Hil ■ 
Tau leptons are identified by their hadronic decays, Th- 
We imposed all the observed leptons e, fj,, m to have 
a transverse momentum larger than 15 GeV, and to be 
within the acceptance of the ATLAS and CMS detectors 
(|ry| < 2.4). The search region to detect experimentally 



SSD events is defined by two additional variables: the 
transverse hadronic energy Ht = J2jet PTjet > determined 
by all the quarks and gluons with px > 40 GeV within 
the detector acceptance and E™ lss , the missing trans- 
verse energy determined by all the undetectable particles 
(neutrinos and LSP) in the event and the visible particles 
outside the detector acceptance. On top of the lepton se- 
lection we require H T > 450 GeV and E™ iss > 120 GeV. 
The SSD detection efficiency is simulated using the for- 
mulas given in [H, [l|| , where the probability for an event 
to pass the selection cuts is given as a function of the Ht, 
i?™ lss and the pr of the two same sign leptons. Events 
with three or more selected leptons arc discarded. Some 
of the SSD events may even originate from a multi-lcpton 
final state, typically with 3 or 4 leptons, where only two 
of them satisfy the selection criteria. The introduction 
of the kinematic cuts (especially those on pr) lessen the 
accuracy of the approximations in Eqs. ([3]) (j4|), which 
have however been verified to hold when no cuts are ap- 
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Figure 3. Curves in the (w.1/2, tan/3) plane showing the luminosity in fb _1 (boxed values) necessary to measure the observables 
iVi (left panel) and N2 (right panel) with ± 50% at the y/s = 14 TeV run, given the CMS acceptance cuts and efficiency model 
discussed in the text. 



plied. We emphasize that the results of the simulation 
presented here are obtained considering all the produc- 
tion mechanisms and without imposing any selection on 
the decay chain that lead to SSD. 

We thus build the contour maps of Ni and N2 in the 
( TO i/2i tan 0) plane that are shown in left panels of Fig.[3J 
yfs = 8 TeV, and right panels, yfi = 14 TeV. Fig.[3Jshows 
the luminosity necessary to achieve an accuracy of 50% 
in the measurement of the observables Ni (left panel) of 
N2 (right panel) at y/s = 14 TeV. We consider here only 
the statistical error. Systematic errors can be neglected 
because N\ t 2 are ratios of observed yields measured using 
the same selection criteria and obtained in the same ex- 
perimental conditions. With the y/s = 14 TeV run that 
will follow after the LHC upgrade, with the possibility 
of accumulating up to 100 - 1000 fb _1 of luminosity, will 
offer the possibility to probe large portions of the pa- 
rameter space in the plane (7711/2, tan/3) up to values of 
m i/2 ~ 900 GeV. We do not show the corresponding of 
Fig. [3] for y/s = 8 TeV because the production cross sec- 
tions are smaller and the planned 20 fb _1 are not enough 
to get N±,2 with the same accuracy. The rather small dif- 
ferences that can be seen in Fig. [2] in the values that iVi^ 
take at the two center of mass energies and in Fig. [3] in 
luminosity curves, confirm our expectation that to a first 
approximation the observables iVi 2 are independent of 
the center of mass energy of the collider. 

Comments and summary. Both the CMS and AT- 
LAS collaborations have released results that hint for a 
possible evidence of a Higgs with mass around 124-126 
GeV in the first ~ 5 fb _1 of data obtained with y/s = 7 
TeV (Him. The light Higgs mass in the MSSM receives 
a large contribution from radiative corrections thus rep- 
resents a crucial quantity to test any SUSY model. Up- 
dated analysis of the CMSSM parameter space including 
the new Higgs data [23[ show that Aq /= 0, large tan/3 
and heavy SUSY spectrum are now generally preferred. 



In the case of the confirmation of the discovery it would 
be interesting to make detailed study of the SSD yield 
expected from the smaller parameter space, compared to 
the one explored here, compatible with that Higgs mass. 

In our analysis we have excluded those model points 
found in the right edge of the strips at high m\/2 in which 
the mass difference between the f\ and the LSP (x?) is 
less then the mass of the tau lepton (m T w 1.7 GeV). In 
this case the two-body decay fi — > Xi T is forbidden. The 
?i decays only into suppressed three body final states, 
and is a long lived charged particle that have been pro- 
posed as a solution to the lithium problem I24|. Once 
produced they can decay outside the detector [25|]. The 
sudden reduction of the number of tau SSD N(tt) is re- 
flected in much higher numerical values of N\ and N2 (in 
particular those of A2). For a recent study connecting 
long-lived staus and SSD in models with the gravitino as 
the dark matter candidate see Ref . [26[ . 

In summary, we have shown that if an excess of SSD 
relative to the SM yield is observed then a measurement 
of the proposed observables A12 within a given accuracy 
will allow to pin down a given portion of the stau co- 
annihilation region (tor) of the CMSSM parameter space 
in the (m^, tan/3) plane. This will in turn also give 
access to other informations. For example if we restrict 
the possible values of miw to a certain range then this 
will of course also restrict the possible values of mo since 
these two parameters are related according to the WMAP 
strips of Fig. [TJ The proposed method can evidently be 
extended to all SUSY models predicting events with final 
states containing SSD or to different extensions of the 
standard model that contain Majorana particles such as 
models with weak scale heavy Majorana neutrinos [27j . 
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Note added. After submitting this paper for publi- 
cation both the ATLAS and CMS experiments have re- 
ported new evidence at 5er level for a scalar particle of 
mass around 125 GeV compatible with the Higgs bo- 
son [28|- As already discussed above, and in view of 
the newly reported experimental evidence on the Higgs 



mass, the choice of the trilinear scalar parameter Aq = 
is not allowed any longer since the radiative corrections 
needed to achieve = 125 GeV are driven by the cou- 
plings of the third generation sfermions and especially so 
the Higgs-stop-stop trilinear coupling A t must be large, 
implying large and negative values for Aq. The allowed 
parameter space of the CMSSM in light of the Higgs dis- 
covery is being still more severely constrained and ac- 
cording to some authors this model is already disfavored 
by the data, while others authors do not agree with this 
conclusion. We refer the reader to Refs. [23[ [29[ for a list 
of recent studies showing these different points of view. 
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